Kondo effect and conductance of nanocontacts with magnetic impurities 
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We study the impact of dynamical correlations on the electronic structure and coherent transport 
properties of Cu nanocontacts hosting a single magnetic impurity (Ni,Co,Fe) in the contact region. 
The strong dynamical correlations of the impurity 3d-electrons are fully taken into account by 
combining density functional calculations with a dynamical treatment of the impurity 3d-shell in 
the one-crossing approximation. We find that dynamical correlations give rise to the Kondo effect 
and lead to Fano features in the coherent transport characteristics similar to those observed in 
related experiments. 
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The development of nanoscale spintronics devices 
based on single molecules and atomic size junctions con- 
taining magnetic atoms is a fascinating and challeng- 
ing field of research at the moment An important 
contribution to the electronic structure and transport 
properties of these devices comes from the strongly in- 
teracting d- or /-electrons of the magnetic atoms. The 
strong interactions result in dynamical correlations that 
give rise to interesting effects like e.g. the Kondo ef- 
fect. For example, Fano lineshapes [2| observed in the 
conductance characteristics of scanning tunneling micro- 
scope (STM) experiments with magnetic adatoms and 
molecular complexes on metal surfaces 0, 3] are the re- 
sult of Kondo resonances at the Fermi level [H, @] . Re- 
cently, Fano lineshapes have also been observed in the 
non-linear conductance characteristics of chemically ho- 
mogeneous nanocontacts 0] made from ferromagnetic 
transition metals (Ni, Co, Fe) [||. 

State of the art calculations of the conductance and 
current through atomic- and molecular-size conductors 
consist in combining ab initio electronic structure calcu- 
lations on the level of density-functional theory (DFT) 
with the non-equilibrium Green's function (NEGF) tech- 
nique [§| . This methodology works quite well for metallic 
nanocontacts Q predicting zero-bias conductances that 
are in general in good agreement with experiments 11|. 
However, static mean-field methods like DFT cannot de- 
scribe dynamical electron correlations. Thus the DFT 
based ab initio transport methodology is not capable 
of describing the Fano-Kondo lineshapes observed in 
STM studies of magnetic adatoms on surfaces 

In order to explore the impact of strong dynamic 
correlations on the transport properties of atomic- and 
molecular-size conductors, we study Cu nanocontacts 
hosting magnetic impurities in the contact region. Such 
a system could also be realized experimentally with e.g. 
the break junction technique using alloys containing 
magnetic atoms like e.g. Cupronickcl. To this end we ex- 
tend the established DFT based ab initio quantum trans- 
port methodology to incorporate dynamic electron cor- 
relations by adapting the LDA+DMFT method HH to 



the case of a single magnetic impurity in a nanocontact. 
While the strong dynamic correlations of the impurity 
d-electrons are fully taken into account, the rest of the 
system is described on a static mean-field level in the lo- 
cal density approximation (LDA) to DFT. Other recent 
approaches to include dynamic electron correlations in 
the ab initio description of quantum transport are based 
on the GW approximation (GWA) [3| or the three-body 
scattering formalism (3BS) [HI]. While the GWA is only 
suitable for weakly correlated systems due to the pertur- 
bative treatment of the electron-electron interactions, the 
3BS is in principle capable of describing more strongly 
correlated systems as it goes beyond perturbation the- 
ory. However, the 3BS does not provide a satisfactory 
solution of the Anderson impurity problem since the lo- 
cal correlations are not taken into account properly. In 
contrast, in our method both the strong Coulomb interac- 
tions between the impurity <i-electrons and the resulting 
local correlations are taken into account properly. 

We consider a single magnetic impurity bridging the 
tips of two semi-infinite Cu nanowires of finite width 
grown in the (001) direction as shown in Fig. Q] We 
divide the system into three parts as shown in the right 
panel of Fig. [TJ Two semi-infnite leads L and R, and 
the central device region (D) which contains the central 
magnetic impurity with the strongly interacting 3d-shell 
(d), and the tips of the two electrodes. The device also 
contains a sufficient part of the semi-infinite leads so that 
the two leads L and R are sufficiently far away from the 
scattering region and the electronic structure of the leads 
has relaxed to that of bulk (i.e. infinite) nanowires. The 
effective one-body Hamiltonians of the device region and 
leads are obtained from DFT calculations on the level 
of LDA. Here we use the supercell approach [l6| to ob- 
tain the effective Kohn-Sham (KS) Hamiltonians of each 
part of the system prior to the dynamical treatment of 
the impurity d-shell and the transport calculations. The 
electronic structure of the device region is calculated with 
the CRYSTAL06 ab initio electronic structure program 
for periodic systems [TtJ by definining a one-dimensional 
periodic system consisting of the device region as the 
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FIG. 1: (Color online) Left: Atomic model of a Cu nanocon- 
tact with a magnetic impurity (red) in the contact region. 
Right: Division of system into left (L) and right (R) elec- 
trode, and central device region (D) containing the magnetic 
impurity (red) hosting the strongly correlated d-orbitals. 



unit cell. The device Hamiltonian H^i is then obtained 
from the converged KS Hamiltonian of the unit cell of the 
periodic system. In the same way, the unit cell Hamilto- 
nians H^ R and hoppings V L / R between unit cells of the 
left and right leads can be extracted from calulcations of 
infinite nanowires with finite width since the electronic 
structure in the semi-infinite leads has relaxed to that of 
an infinite nanowirc. In the LDA calculations we employ 
a minimal basis set plus effective core pseudo-potential 
that takes into account only the 4s, Ap and 3d valence 
shells of the Cu atoms and the magnetic impurity [la |. 

The strong electron correlations in the 3d-shell of 
the magnetic impurity are captured by adding a 
Hubbard-like interaction term to the one-body Hamil- 
tonian within the correlated subspace d: jixj = 
\HUi ] uc\ <J1 c { J<J2 ci C!2 c kai (Einstein sum convention). 
Uijki are the matrix elements of the effective Coulomb 
interaction of the 3d-electrons which is smaller than the 
bare Coulomb interaction due to the screening by the 
conduction electrons. In the spherical approximation all 
matrix elements Uijki can be calculated from the Slater 
integrals F°, F 2 , and F 4 which are related to the aver- 
age Coulomb repulsion U between electrons and to the 
Hund's rule coupling J by F = U, F 2 = (14/1.625)J, 
and F A — 0.625F 4 [13j. For 3d transition metal elements 
in bulk materials the repulsion U is around 2-3 eV and 
J is around 1 eV [19j. Due to the lower coordination 
of the contact atoms the screening of the direct interac- 
tion is reduced compared to its bulk value. Here we take 
U = 5 eV and J = 1 eV, but we have checked that the 
results do not change much when U is varied between 4 
and 6 eV. 

The Coulomb interaction within the correlated 3d sub- 
space has already been taken into account on a static 
mean- field level in the effective KS Hamiltonian of the de- 
vice. Therefore the KS Hamiltonian within the correlated 
subspace Hd has to be corrected by a double-counting 
correction term, i.e. Hd = H^ s — Hd c - Here we use the 
standard expression, Hd c = [U(Nd — |) — \ J(Nd~ 1)] xl<j 
where Id is the identity matrix in the d subspace, and Nd 
is the occupation of the impurity 3d-shell [13] . 

The central quantity is the Green's function (GF) of 



the device region: 

G D = (cj + ^ - H D + H dc - £ d - £ L - Er)- 1 (1) 

where [i is the chemical potential. El and Er are self- 
energies that describe the coupling of the device to the 
semi-infinite leads L and R, respectively. These can 
be calculated from the effective one-body Hamiltonians 
of the leads by iteratively solving the Dyson equation 
S l /r = V L/ R(cj+/i-H L/R -S L/ R)- 1 v[ /R . E rf is the lo- 
cal electronic self-energy that describes the dynamic elec- 
tron correlations of the impurity 3d-electrons. In order 
to calculate Ed, the generalized Anderson impurity prob- 
lem given by the impurity 3d-shell has to be solved. The 
impurity problem is described by the projection Pd of the 
GF (Q} onto the correlated subspace d: Gd = PdGnPd 
which can be written as 



G d (w) = (w + n - H rf - E«,(w) - A d (w))" 



(2) 



where we have introduced the so-called hybridization 
function Ad which describes the hybridization of the 
impurity electrons with the conduction electrons. The 
hybridization function can be calculated from the pro- 
jection Pd of the uncorrelated device GF G° D = 
(lj + fJ, — Hd + Hdc — El — S^) _1 onto the correlated 
subspace d [16|, i.e. from G° = PdG° D Pd- Solving eq. 
© for A d and using that [Gd]" 1 = [G^]- 1 - E d we 
obtain: 



A d (w) 



(3) 



The hybridization function Ad [16[, the Coulomb re- 
pulsion U, the Hund's rule coupling J, and the impurity 
levels ed,i — (Hd)ii are the relevant parameters for solv- 
ing the impurity problem. Here we employ the so-called 
One-Crossing - Ap proximation (OCA) to solve the impu- 
rity problem [13l . [2fj which is particularly well suited for 
the Kondo regime. 

The current through a strongly interacting impurity 
can be calculated exactly by the Meir-Wingreen formula 
2l| . However, for low temperatures and small bias volt- 
ages this expression is well approximated by the much 
simpler Landauer formula j22|: I(V) — ^ X J Q e dujT(Lu) 
where T(oj) is the Landauer transmission function and 
where we have assumed an asymmetric voltage drop V 
about the device region j23j. Thus the conductance 
is simply given by the Landauer transmission function: 
Q{V) = J£(V) = 2|! x T(eV). The latter can be 
calculated from the device Green's function: T(uS) = 
Tr[T L (u})G\j(u)T R (Lj)G D (u})} where T L/R are the so- 
called coupling matrices which describe the coupling to 
the leads, and can be calculated from the lead self- 
energies by T L/R = i(S L / R - E^ /fl ). 

Fig. [2] shows the result of our LDA+OCA calcula- 
tions for three different impurities: Ni, Co, and Fe, for 
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FIG. 2: (Color online) PDOS of d-orbitals calculated with the 
LDA+OCA method for different magnetic impurities in Cu 
nanocontact for the geometry shown in Fig. [JJ (a)-(c) PDOS 
near Fermi level for a Ni (a), Co (b) and Fe (c) impurity 
at different temperatures, (d) Comparison of the PDOS of 
the three impurities on a larger energy scale at T = f 20 K. 
(7 = 5 eV and J — 1 eV in all cases. 
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FIG. 3: (Color online) Conductance calculated with 
LDA+OCA method for different magnetic impurities in Cu 
nanocontact for the geometry shown in Fig. [JJ (a)-(c) Con- 
ductance Q vs. bias voltage V for different temperatures and 
for small bias, (d) Comparison of LDA+OCA (solid lines) and 
LSDA (dashed lines) conductances for different impurities at 
T = 120 K. U = 5 eV and J = I eV in all cases. 



different temperatures. In all three cases the partial den- 
sity of states (PDOS) for the impurity 3d-electrons shows 
resonances near the Fermi level which are temperature- 
dependent. More precisely, the resonances vanish with 
increasing temperature. This is characteristic for the 
Kondo effect which is usually observed only at low tem- 
perature. The resonances originate from different d- 
orbitals in each indicated by the labels in the 

figures. In the case of Ni, the resonance originates from 
the d^y-orbital. In the case of Co, there are two distinct 
peaks corresponding to two different sets of orbitals: The 
peak that is farther from (closer to) the Fermi level orig- 
inates from the d x 2_ y 2 (d xz ,d yz ) orbital(s). The doubly- 
degenerate d xz , d yz -orbitals are also responsible for the 
resonance in the case of Fe. As can be seen from Figs. 
[3](a)-(c) the corresponding conductances all show Fano- 
like features. Interestingly, in the case of Co, the d x i_ y i- 
resonance does not lead to a corresponding feature in 
the conductance characteristics. This can be understood 
by the so-called orbital blocking: The electron transport 
through certain orbitals can be inhibited by the geome- 
try or symmetry of atomic-size conductors in spite of the 
orbital having spectral weight near the Fermi energy [Io| . 

Tab. [J shows the orbital occupations and effective en- 
ergy levels Id = £d + Re Ad(0) of the impurity 3d-shell 
for different impurity atoms. One can see, that the ef- 
fective energy levels id roughly correlate with the orbital 
occupations except in the case of the d x 2_ y 2 -orbital in Ni 
and Fe. To fully understand the orbital occupations, the 
imaginary part of the hybridization function has to be 
taken into account. The imaginary part of describes 
the broadening of the d-orbitals due to the coupling to the 
conduction electrons. It turns out that the d x i_ y i -orbital 



is by far the most localized orbital of all the d-orbitals 
i.e. the broadening is very small compared to the other 
d-levels This explains why the d 2; 2_ l ,2-orbital in Ni 
and Fe is only half-filled despite its low effective energy 

id- 

Furthermore, we see that for Ni the d^-orbital which 
gives rise to the resonance near the Fermi energy is almost 
completely filled. Hence it does not carry a spin-1/2 and 
therefore the system is not in the Kondo regime, but 
is in the so-called empty orbital regime 24J. A broad 
quasiparticle (QP) peak quite close to the Fermi level 
appears which — as in the Kondo regime — arises from 
the Fermi liquid behaviour at low temperatures. The 
temperature dependence of the QP peak is qualitatively 
similar to the Kondo regime, i.e. the peak broadens with 
increasing temperature and disappears above a critical 
temperature, called the Kondo temperature Tk. For the 
Co and the Fe impurity, the doubly-degenerate orbitals 
d xz and d yz that give rise to the QP resonance at the 
Fermi level are occupied by three electrons (filling 3/4) 
and thus carry a spin-1/2. Therefore in the case of Co 
and Fe, the system really is in the Kondo regime. 

We have estimated the Kondo temperatures Tk dis- 
played in Tab. U from the width of the QP peaks at 
low temperature. The Kondo temperatures follow the 
same trend, namely T K (Ni) > Tk(Co) > Tk(Fe), which 
is also observed in STM experiments with adatoms on 
metal surfaces |25j and in pure transition metal nanocon- 
tacts Moreover, the Kondo temperature for Co agrees 
quite well with the Tk estimated from recent STM ex- 
periments with Co adatoms on Cu surfaces in the con- 
tact regime Q ■ The high Kondo temperatures of about 
200 K for Co and Fe imply a strong antiferromagnetic 
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N d 


8.42 


7.04 


6.02 


Tk [K] 


~ 500 


~ 225 


~ 160 



TABLE I: Orbital occupations n d , effective energy levels 
£d = e d + Re Ad(0) of impurity d- levels relative to d 3z 2_ r 2- 
level, total occupation N d of the impurity 3d-shell and Kondo 
temperature Tk estimated from the width of the resonance 
closest to the Fermi level for each of the three impurities. 
U = 5 eV, J = 1 eV. 

coupling J sc i between the conduction electrons and the 
impurity <i-electrons giving rise to the Kondo effect since 
Tk oc exp(— 1/ J s dPo) where on is the conduction elec- 
tron DOS at the Fermi level [24J. This strong antiferro- 
magnetic coupling might explain why the Kondo effect is 
observed in ferromagnetic nanocontacts despite the fer- 
romagnetic coupling to the bulk electrodes [8j. 

Finally, in Fig. [3)Jd) we compare the results obtained 
with the LDA+OCA method at low temperature with re- 
sults obtained from DFT calculations on the level of the 
local spin density approximation (LSD A). For Ni (red 
lines) the effect of including dynamic correlations is only 
moderate. Thus for Ni the static mean-field description 
given by LSDA is a reasonable approximation to the fully 
correlated description by the LDA+OCA method, but 
at the cost of breaking the spin symmetry. This can 
be understood by recognizing that LSDA usually gives 
reasonable spectra for the empty orbital and mixed va- 
lence regimes. In contrast for Co (green lines) and Fe 
(blue lines) taking into account dynamical correlations 
changes the conductance substantially. The dynamic 
correlations open a gap in the 3d-shell thereby taking 
away spectral weight from the Fermi level, leaving only 
the Kondo resonance (with small spectral weight) at low 
temperatures. Consequently, the conductances predicted 
by LDA+OCA are considerably lower than the conduc- 
tances predicted by LDA and LSDA. 

In conclusion, we have extended the established DFT 
based ab initio transport methodology for nanoscopic 
conductors to include dynamic electron correlations. We 
find that nanocontacts hosting a magnetic impurity show 
strong dynamical correlations which give rise to quasi- 
particle resonances at the Fermi level and corresponding 
Fano features in the conductance- voltage characteristics. 
Our findings agree well with experiments measuring the 
conductance through Co adatoms on metal surfaces in 
the contact regime. Moreover, our results shed some 
light on the recent observation of the Kondo effect in 
ferromagnetic nanocontacts. 
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